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Abstract. I present a short overview of the NLO QCD calculations available for 



^ ■ deep-inelastic production of heavy quarks. 
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1. Introduction 

Charm quarks produced in deep-inelastic scattering (DIS) have been identified in sizable 
numbers now by the HI [1] and ZEUS [2] collaborations at HERA, and considerably 
more charm (and bottom) data are anticipated. At the theoretical level the reaction has 
already been studied extensively. In the framework where the heavy quark is not treated 
■ as a parton, leading order (LO) [3, 4] and next-to-leading order (NLO) [5, 6] calculations 
of the inclusive structure functions exist. Moreover, LO (AROMA, RAPGAP) [7, 8] 
and NLO (HVQDIS) [9, 10, 11] Monte-Carlo programs, allowing a much larger class of 
observables to be compared with data, have been constructed in recent years. Overall, 
the NLO QCD description agrees quite well with the HERA data. Here I shall give a 
very brief overview of these NLO calculations. 

Charm quarks are produced in DIS via the reaction 

^ ! e ± (p e ) + P(p) -> e ± (p e -q) + X[Q, Q] , (1) 

where P(p) is a proton with momentum p, Q is a heavy quark with momentum p\ 
(jp\ = m 2 ) and X is any hadronic state allowed, containing the heavy quark-antiquark 
pair. Its differential cross section may be expressed in general as 



d 2+n a 2ira 



dxdQ 2 Hi dVi x Q 4 
where 



rp Jn rp 



(2) 



Q = -q , x=- , y = . 3 

2p ■ q p-p e 

The Vi stand for kinematic variables related to the heavy quarks. Examples are the 
transverse momentum of the heavy quark, the rapidity difference between the heavy 
quarks, etc. 
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2. Inclusive and single-charm inclusive production 

The least difficult cross section to measure is the heavy quark inclusive cross section, 
expressed via (2) in the inclusive structure functions F 2 and F L . These were calculated 
to NLO in [5]. The results are parametrized as 

Q 2 a s /■zmax dz 



+ ^/rf[4/ s (p^)(4>C.n^) 



-nm? 



+ E [4 t^Ml + eg In + el, fA ^)dgl ] ] , (4) 



where k = 2,L and the upper boundary on the integration is given by z max = 
Q 2 /(Q 2 + 4m 2 ). The parton densities f\(x, fi 2 ) , (i = g,q,q) are explicitly identified. 
The scale fi is the mass factorization and renormalization scale. The MS coefficient 

functions c^-fa, , cj^fa , (j = 9,Q,Q',l = 0,1) and djJifaO* (* = 9> 95^ = °> X ) 
depend on the scaling variables rj — (s — Am 2 ) /Am 2 and £ = Q 2 /m 2 , with s the square 
of the cm. energy of the virtual photon-parton subprocess. This implies that in (4) 
z = Q 2 /(Q 2 + s). 

In (4) we distinguish the coefficient functions with by their origin. The 
c k!i(Vi O'Cfc'k 7 ?) originate from those partonic subprocesses where the virtual photon 
is coupled to the heavy quark, whereas the d^^r],^) correspond to the subprocess where 
the virtual photon interacts with the light quark. Thus the former are multiplied by the 
charge squared of the heavy quark e 2 H , and the latter by the charge squared of the light 
quark e\ respectively (both in units of e). Only the terms proportial to e 2 H contain the 
gluon density. 

To obtain numerical results for the inclusive cross section, it is better to use instead 
of the original, rather long expressions [5], the much faster parametrized form [6]. The 
lowest order term contains only the gluon density. Light quark densities only come in 
at next order, contributing only about 5%. This is the reason F 2 (x,Q 2 ,m 2 ) is used 
in global analyses to constrain the gluon density. Besides the gluon density, the main 
source of theoretical uncertainty in F 2 (x, Q 2 , m 2 ) is the value of the charm quark (pole) 
mass, rather than the scale /i. 

The calculation that lead to (4) also yielded the single heavy quark differential cross 
section, with V\ = and V 2 = y Q [12]. These distributions are best generated using 
the HVQDIS program, described in the next section. 

NNLO estimates based on soft gluon resummation are given in Ref. [13], for 
inclusive and single heavy quark inclusive cross section. 



3. Fully differential charm production 



A NLO calculation also exists for the fully differential cross section in (2) [14]. 
Maintaining full different iality required a complete recalculation of the matrix elements, 
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carefully eliminating intermediate divergences via the so-called subtraction method. The 
results are encoded in the program HVQDIS [10] \. The program can compute, to NLO, 
experimentally visible cross sections, which are in principle better for comparison with 
theory than fully inclusive ones. It returns parton kinematic configurations and their 
corresponding weights, accurate to 0(aa 2 s ). The user is free to histogram any set of 
infrared-safe observables and apply cuts, all in a single histogramming subroutine. 

Additionally, one may study heavy hadrons using the Peterson et al. fragmentation 
model. Detailed physics results from this program are given in [11]. 

HVQDIS has been used extensively in experimental analyses. As is shown elsewhere 
[16] in these proceedings, it reproduces the data very well indeed, except for L>*'s at low 
Pt and large pseudorapidity, where there are more events than HVQDIS would predict. 
This is possibly due to remnant beam drag effects distorting the pseudorapidity spectrum 
to larger values. For the case of charm photoproduction this was investigated in [17]. 

A more extensive overview of the NLO calculations and the phenomenology of DIS 
charm production can be found in [18]. 
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f See [15] for a description of the most recent changes. 



